Statistical properties of localisation— delocalisation transition in one dimension 
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We study a one-dimensional model of disordered electrons 
(also relevant for random spin chains), which exhibits a de- 
localisation transition at half-filling. Exact probability distri- 
bution functions for the Wigner time and transmission coef- 
ficient are calculated. We identify and distinguish those fea- 
tures of probability densities that are due to rare, trapping 
configurations of the random potential from those which are 
due to the proximity to the delocalisation transition. 

PACS No: 



The Anderson transition in dimensions D < 3 has re- 
cently attracted a renewed interest in relation to such 
systems as random antiferromagnetic spin chains ^ and 
high mobility Si MOSFET's g. Experiments on both 
systems could not be accounted for in the standard scal- 
ing theory of localisation j|] . 

The simplest disordered model known to exhibit metal- 
lic behaviour, in spite of being one-dimensional, is the 
random-hopping model, 



(1) 



where i„ > are random variables, with n-independent 
average, t, and c„ annihilates a spinless fermion 

at site n. This ID model has a single delocalised state 
at the middle of the band, e = 0. This is an interesting 
example of the failure of the scaling theory of localisation. 
Moreover, this model has many common features with a 
wide class of random spin chains; such as the spin-1/2 
random Heisenberg chain H ~ JnSn ■ Sn+i, where 
Jn > are randomly distributed. (The XX version of 
the latter model is, in fact, exactly equivalent to (|^), 
upon the Jordan- Wigner transformation.) 

For the random hopping model, a great deal is known 
about such self-averaging quantities as the total density 
of states. More recently, some of the correlation func- 
tions have also been calculated. However, the behaviour 
of probability distributions in the proximity of the de- 
localisation transition is virtually unexplored. It is the 
purpose of this letter to address this issue. 

In the continuum limit, model (^) becomes what is 
known as the random-mass Dirac model (see e.g. [||J^): 

n = -i[R^d^R ~ L^d^L ]~im{x)[R^L -L^R] (2) 

where R and L are the chiral components of the electron 



field operator. The derivation of the continuum limit as- 
sumes weak disorder such that t„ = t + 5tn (the Fermi ve- 
locity, vp, associated with t is set to 1). It is the staggered 
component of the random hopping, Stn — *■ {—l)""m{x), 
which enters into the continuum theory. In field-theoretic 
language, this corresponds to a random mass. 

It was found by Dyson in 1953 g, that the average 
electron density of states for a model equivalent to (|l|) di- 
verges at the middle of the band: p(e) ^ l/(e| Inep). By 
the Thouless relation 0| , such a density of states implies 
a divergent localisation length Ae ~ | lne|. The criticality 
of the model at half-filling was established by Gogolin 
and Mel'nikov In particular, they found that model 
(^) has a finite conductivity in contrast with the Mott 
law in the standard localised regime. Calculations of Ref. 
Q first indicated that it is not the Thouless length A^, 
but rather the length Ig ~ In^ e, which is likely to gov- 
ern the correlation functions. The role of the length Z,. 
was later clarified by means of the real-space renormal- 
isation group method by Fisher The cross over for 
the electron Green function between the spatial regions 
a; <C ?e (critical regime) and <C a; (localised regime) was 
discussed in recent publications [Q and Q by means of 
the Berezinskii technique and the Efetov supersymmetry 
method, respectively. 

It seems natural to undertake the next step and in- 
vestigate the probability distribution functions for the 
above systems. To this end, we adopt a simple math- 
ematical technique, based on the recent observation by 
Shelton and Tsvelik that, at e = 0, the random mass 
Dirac model can be formulated as a one-particle quantum 
mechanical problem. Indeed, if we introduce the combi- 
nations x± = (RT L)/^/2 of the chiral components of 
the electron field, then the Dirac equation becomes 

[dx T m{x)]x± = iexT ■ 

This equation decouples at e = 0, thus admitting the 
(unnormalised) solutions of the form 

X±ix) ^ e±^(-). Vix) = r dym{y) (3) 

Jo 

As it is customary in the literature, we assume that 
the random mass m[x) is i5-correlated in the real space 
and Gaussian: P„ig[m{x)\ ^ exp[— (m(a;) — mg)^ /2(7]. 
Here g characterises the disorder strength and the mass 
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has a mean value, mo- Thus, solution (^) for the wave- 
functions is nothing but an exponential of a random walk 
with a drift term, mQX. 

Wave-function (||) allows us to investigate the distri- 
bution of several physical quantities in a relatively sim- 
ple manner. We start by considering the so-called re- 
laxation time, which was introduced in the context of 
the scattering theory by Wigner [Q, hence also known 
as the Wigner time. Physically this is the time spent 
by a wave packet inside a scattering region, and it can 
be formally defined as the momentum derivative of the 
scattering phase shift. Let the effect of the disorder be 
confined to the segment [0, L], which is the scattering re- 
gion in our problem. For the sake of concreteness, we 
assume that the electrons can not leave the sample on 
the left-hand-side (x < 0), so that they are scattered off 
the segment [0, L] on the right. We thus impose a van- 
ishing boundary condition on the left, X-(0) = (this 
condition corresponds to suppressing a site in the lattice 
formulation, (|l|)). The boundary condition on the right 
is simply the continuity of the Dirac wave-function at 
X = L. In this case, the calculation of the Wigner time 
involves the scattering phase picked up by the incident 
left-moving wave in the process of scattering from the 
disordered segment, so that (x > L): 



T, = d0/de, e{e) = -Un [e-'^'"'' R{x) / L{x)\ 



(4) 



We are mainly interested in the relaxation time at half- 
filling, i.e. tq. In order to find this, we perturb the Dirac 
equation in e around the e = solution (^. Perform- 
ing this simple calculation we found the following exact 
expression for the Wigner time as a functional of the dis- 
order [d 



(5) 



In accordance with |l2| , the relaxation time is related to 
the total charge, Qo = /g^ '^2;(|i?oP + |ioP)j and therefore 
To — Qo — 2L. 

As such, expression (^) does not supply much infor- 
mation. Indeed, the physical content of the problem is 
revealed by the probability density of the relaxation time 
^ko] — ^ Ti3\V])) . We observe that the Laplace 

transform of this density, P[A] = droe~'^'^°7^[To], can 
be represented as a path integral 



P[A]~ / W(x)F„,jV^]e-^^«[^l 



(6) 



In the above formula, the integration is taken along all 
paths starting at V{x — 0) and ending at V{x = i) = Vq. 
By construction, see (^), the starting point is V^(0) — 0, 
while the ending point Vq is arbitrary, so the path integral 
in (|^) involves an additional integration over all possible 
end points Vq. Fixing the normalisation, P[0] = 1, and 



performing a convenient shift of the integration variables, 
we ultimately obtain 



P[\] 



dVoKx{Vo,0;L) 



(7) 



(For the time being we have set mo — 0, so as to be 
at the criticality.) Here K\{Vi,V2; x) is the imaginary- 
time propagator for a quantum mechanical particle char- 
acterised by the action 



SxlV] = 



dx 
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(8) 



The system described by this action is known as Liou- 
ville quantum mechanics [|l4| . The Schrodinger equation. 



corresponding to the action 
malised) solution jlj] 



), has the following (nor- 



(9) 



with the energy — 2A + {Ki-y is the MacDonald 

function). Constructing the propagator in a standard 
manner by making use of the complete set (^, integrating 
over Vq, and performing the inverse Laplace transform, 
we obtain the exact Wigner time probability density 



nro] 



2l/2g7rV8gL 



dt cosh t cos 



exp[-cosh^t/{gQo) ^t^/{2gL)] . (10) 



It is convenient to represent the result in terms of the 
positively definite quantity (total charge) Qq — tq + 2L. 
According to the general theory of Ref. the delay 
time can be negative (for repulsive potentials) but there 
is a lower bound, which is —2L in our case. 

We note the following interesting limiting cases of the 
formula 

• Consider first the probability of long time delays, 
Tq 00 (and therefore Qq — > 00). The i-integral in 
( pO| ) is convergent for all Qq. It is therefore tempting 
to expand the exponential in the integrand in powers of 
l/Qo- It is easy to check, however, that all the coefficients 
of such an expansion identically vanish. It follows that 
( pO| ) has an essential singularity at tq = cx). The nature 
of this singularity can be determined by first neglecting 
the l/Qo term in the exponential and then simulating 
its effect by cutting off the i-integration at large times 
to (1/2) lii((7ro), when the term becomes of the order 
of unity. Then, within the leading logarithmic accuracy, 
we obtain: 



PN-exp [-\n'{gTo)/8gL] 



(11) 



This formula is valid when the factor in the exponential 
is large, i.e. In^ tq ^ L. As a simple application to 
random spin chains, we consider long-time relaxation of 
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the magnetisation M{t) inside a finite segment of length 
L. Due to (|ll|) we find that M{t) - exp [-ln(t)/8gL]: 
a very slow decay. 

• Since the applicability of (|l^) involves the system size 
L, the probability distribution for large samples ought to 
be different. Indeed, for a fixed Qo, we obtain: 



1 



-1/9Q0 



^/TrgLQo 



(12) 



The next quantity of interest is the transmission coef- 
ficient at e = 0, which is also proportional to the Lan- 
dauer conductance at half-filling. In order to have a finite 
transmission coefficient T, let us open our sample on the 
left-hand-side, so that the electrons can now leave it at 
X = (where the boundary condition thus becomes the 
same as at a; ^ L). Upon matching the wave- functions 
in the usual way, one finds 



T[V] = l/coshV(L) . 



(13) 



This is a simple formula: the transmission coefficient in 
not a functional of the entire random walk trajectory 
(as the Wigner time is) but only a function of the end 
point. Therefore the probability distribution P[T] can be 
found in an elementary way without using the Liouville 
mechanics. We obtain 



P\T\ = - — , exp 



1 

25! 



cosh"^(l/\/T)j I (14) 



where < T < 1. This is an intriguing distribution 
function, plotted in Fig.l. It has the following properties: 

• The probability of a small transmission is: 

P[r]^exp[-ln2(l/r)/8.gL] (15) 

• The function P[T] has a low-transmission peak at Tq 

• There is an integrable divergence close to the perfect 
transmission (T — 1): 



P[T] ~ ^8/^5^(1 - T) 
The mean transmission coefficient is given by 



(16) 



(17) 



This formula is asymptotically exact as i — > cxd (we also 
verified this result by an independent calculation using 
the Berezinskii technique [T5|]). This result is surprising. 
Since T ^ cr/L, a being the dc-conductivity, the 1/\/L 
behaviour of the average transmission coefficient suggests 
that the system is even more metallic than what one 
would expect in an hypothetical case of a weakly dis- 
ordered ID metal (hypothetical because any weak dis- 
order is supposed to lead to localisation). It is worth 
noting that the transmission coefficient is not a good 
scaling variable in the limit L —^ 00. So, from Eq.(14), 



one finds that the variance of the transmission coefficient 
= ((T - (T))2) - 1/VX, implying that the width 
of the distribution, normalised to the mean transmission 
5/{T) ~ L*, diverges when L — > 00. On the other hand, 
the logarithm of the transmission coefficient is a good 
scaling variable in the above sense, with average 

(log (^^^) = a^/^, 

a being a positive numerical coefficient. This conclusion 
was already reached in 1980 by Anderson et al. in study- 
ing a generic ID disordered system |lq| . They introduced 



the scaie conductance, Tt 



3{iogiT instead of the av- 
erage conductance. In our model we reach a similar con- 
clusion. However, unlike the standard case, we find that 
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Ttyp = e~"V^^ behaves quite differently from the aver- 
age transmission, and both are different from the peak in 
the probability distribution Tq ~ e^*^^. 

There is a difference between our 1 / result for the 
Landauer conductance and the finite conductivity found 
for an infinite system in Ref. A possible physical ex- 
planation is that this difference is due to resonant scat- 
tering processes, which enhance the probability of near- 
perfect transmission in finite samples but are absent in 
infinite systems (this is characteristic for a critical sys- 
tem; otherwise, the T ^ 1 divergence of P[T] is expo- 
nentially suppressed in the sample length). However, this 
question requires further studies. 

Both calculations for V[to] and P[T] can be generalised 
to the off-critical case when mo is non-zero. The details 
will be given in an extended article JTst . Here we only 
quote the results. Formulae ([l|) and (|l^) are not af- 
fected. The power law in the denominator of (112) changes 



to 1/Q, 



l+mo/g 




Both formulae (|T|) and (|T|) acquire a 
suppressing factor, exponentially small in the parameter 
moL. 
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FIG. 1. Distribution of the transmission coefficient with 
gL = 1.25 and average mass zero (solid line) and moL = 0.5 
(dashed line). 
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Notice that the probabiHty of long time delays, (11), 
as well as the probability of a small transmission, (15), 
follows the so-called log-normal law. To our knowledge, 
the log-normal tails of the distribution functions in one- 
dimensional disordered systems were first obtained in 
Ref. 1 16 1 by means of a scaling argument. Their ex- 
istence was rigorously established by Mel'nikov via the 
Bereziskii technique [|l^. In 2 -I- e dimensions, the log- 
normal tails were found in Ref. ||l^ , which was thought to 
be a signature the Anderson transition. The discussion of 
the log-normal tails was recently revived by Muzykantskii 
and Khmelnitskii |l9| (see also ||2^) who gave a simple 
derivation based on a specific saddle-point approximation 
in the supersymmetric cr-model. (For a recent collection 
of results on the log- normal distributions see |^ .) 

We have shown that the random mass Dirac model 
does posses log-normal tails at the criticality. More- 
over, these tails are unaffected when one moves away 
from the criticality. This is consistent with the inter- 
pretations 1 19 2^ that these tails are due to the so-called 
'anomalously localised' electronic states, which occur in 
rare, trapping disorder configurations. Indeed our result 
( pT| ) can be understood in terms of the 'optimal fluc- 
tuation' concept, discussed in this context in ||l^. The 
optimal fluctuation in our case corresponds to having a 
constant mass mo within the sample. This would ac- 
cumulate charge Qo oc exp(2moi). The probability to 
have such a potential, and therefore such a charge is 
~ exp(— mgL/2(;) ^ exp(— In^ Qo/8gL). 

So, does it follow that the proximity to the delocalisa- 
tion transition (criticality) plays no role? To clarify this 
point we first notice that the log-normal tails are only 
present in a finite system. (These tails do not follow from 
and are not directly related to the so-called multi-critical 
exponents appearing in the wave- functions' statistics for 
the infinite system pjll|].) The Wigner time distribution 
function takes a different form in the limit of a large 
system, (|l2|). This can be interpreted as a 'limiting' 
('equilibrium') distribution in terms of the Fokker-Planck 
equation approach ||l3| and it reveals no trace of the log- 
normal behaviour. It is this distribution which bears the 
signature of the criticality: it ceases to be normalisable 
for Too = thus requiring a long-time cut off. Similarly, 
the log-normal tail appears at low— T in the transmis- 
sion coefficient distribution function, while the criticality 
shows up near the perfect transmission. Indeed, the di- 
vergence of P[T —> 1] is always in place (due to resonant 
scattering processes) but it is exponentially suppressed 
in the sample length, unless the system is critical. As 
a result, at the criticality, the mean transmission coeffi- 
cient is not any more exponentially small but is given by 
the power-law formula (p7|). (Notice, though, that the 
peak of P[T] still exists at the criticality.) It follows that 
the disorder configurations leading to the log-normal tails 
and to the delocalisation phenomenon act independently. 
They affect different domains of the parameter space of 



the problem and show up in the distinct limits for the 
probability distribution functions. 

We are thankful to A.M. Tsvehk, D.E. Khmelnit- 
skii, and B.A. Muzykantskii for interesting conversations. 
M.F. acknowledges support by INFM, Project HTSC. 



[1] M. Fabrizio and R. Melin, Phys. Rev. Lett. 78, 3382 
(1997). For an experimental review on spin-Peierls sys- 
tems see J. P. Boucher and L.P. Regnault, J. Phys. I. 
France 6, 1936 (1996). For doping effects in spin lad- 
ders see M. Tanako, Nature 377, 41 (1995); M. Azuma, 
et al, Phys. Rev. B 55, R8658 (1997). 

[2] S.V. Kravchenko, D. Simonian, M.P. Sarachik, W. Ma- 
son, and J.E. Furneaux, Phys. Rev. Lett. 77, 4938 (1996). 

[3] E. Abrahams, P.W. Andesron, D. Licciardello, and T. 
Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 

[4] M. Steiner, M. Fabrizio, and A.O. Gogohn, Phys. Rev. B 
57, 8290 (1998). 

[5] L. Balents and M.P.A. Fisher, Phys. Rev. B 56, 12 970 
(1997). 

[6] F.J. Dyson, Phys. Rev. 92, 1331 (1953). 

[7] D.J. Thouless, J. Phys. C 5, 77 (1972). 

[8] A.A. Gogolin and V.I. Mel'nikov, Zh. Eksp. Teor. Fiz. 
73, 706 (1977) [Sov. Phys. JETP 46(2), 369 (1977)]. See 
also the review paper A.A. Gogolin, Phys. Rep. 86, 1 
(1982). 

[9] J. P. Bouchaud, A. Comtet, A. George, and P. Le Doussal, 
Ann. Phys. 201, 285 (1990). 
D.S. Fisher, Phys. Rev. B 50, 3799 (1994). 



D.G. Sh elton and A.M. Tsvelik, Report No. :ond-mat/ 



970606C 



[12] E. Wigner, Phys. Rev. 98, 145 (1955). 

[13] A similar expression was found, using a different method, 
in the recent paper A. Comtet and C. Texier, J. Phys. A 
30, 8017 (1997). 

[14] E. D'Hoker and R. Jackiw, Phys. Rev. D 3517 (1982). 

[15] M. Steiner, M. Fabrizio, and A.O. Gogolin (unpublished). 

[16] P.W. Anderson, D.J. Thouless, E. Abrahams, and D.S. 
Fisher, Phys. Rev. B 22, 3519 (1980). 

[17] V.L Mel'nikov, Pis'ma Zh. Eksp. Theor. Fiz. 32, 244 
(1980) [JETP Lett. 32, 226 (1980)]. 

[18] C. Castellani, C. Di Castro, and L. Peliti, J. Phys. A 
19, L1099 (1986); B.L. Affsuler, V.E. Kravtsov, and 
I.V. Lerner, Pis'ma Zh. Eksp. Theor. Fiz. 45, 160 (1987) 
[JETP Lett. 45, 199 (1987)]. See also J.S. Caux, I. Ko- 
gan, A. Lewis, and A.M. Tsvelik, Nucl. Phys. B 489, 469 
(1997). 

[19] B.A. Muzykantskii and D.E. Khmelnitskii, Phys. Rev. B 
51, 5480 (1995); Phys. Rep. 288, 259 (1997) 

[20] V.L Fal'ko and K.B. Efetov, Phys. Rev. B 52, 17 413 
(1995). 

[21] A.D. Mirlin, J. Math. Phys. 34, 1888 (1997). 



4 



